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摘 要 : 月 球 物 理 天 平 动 是 对 月 球 运动 在 空间 摆动 的 描述 。 确 定 月 球 物理 天 平 动 , 可 以 


推测 物理 天 平 动 的 激发 与 耗 散 机 制 , 比如 陨 击 、 月 震 和 核 慢 粘 潇 摩擦 等 , 故 测定 月 球 物理 天 
平 动 对 认识 太阳 系 中 天 体 的 起 源 、 演 化 及 结构 等 , 都 具有 十 分 重要 的 意义 。 本 文 利 用 最 新 的 
INPOP19a 历 表 的 数据 ， 完 成 了 对 历 表 提 取 的 欧 拉 角 到 月 球 物理 天 平 动 的 转换 ， 得 到 的 物 
理 天 平 动 数值 分 别 与 该 系列 历 表 INPOP17a 以 及 被 广泛 使 用 的 DE430 对 比 ， 发 现 不 同 历 
表 的 物理 天 平 动 之 差 存在 着 稳定 周期 。 对 比 历 表 欧 拉 角 的 差别 ， 计 算出 地 心 到 月 面 反射 器 
A15 的 距离 最 大 有 30 厘米 的 差别 ， 此 结果 对 月 球 激光 测 距 的 预报 精度 有 较 大 的 影响 ， 为 后 
续 的 高 精度 测 月 研究 打下 基础 。 研 究 结 果 表 明 ，INPOP19a 最 为 稳定 ， 故 在 月 球 物理 天 平 
动 研究 中 推荐 使 用 INPOP19a。 
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月 球 天 平 动 是 对 月 球 运动 在 空间 摆动 的 描述 ， 它 是 月 球 自 转 ， 公 转轨 道 以 
及 周围 天 体 对 月 球 的 力 的 作用 的 综合 表现 。 基 于 激光 测 月 数据 制作 的 月 球 历 表 
是 进行 月 球 物理 天 平 动 研究 的 数据 基础 ”。 

根据 卡 西 尼 定 则 ， 月 球 的 公转 周期 与 自转 周期 一 致 ， 但 由 于 月 球 质量 的 非 
球面 分 布 ， 月 球 的 自转 并 不 是 均匀 的 。 这 种 与 均匀 旋转 的 偏离 被 称 为 物理 天 平 
动 。 物 理 天 平 动 有 两 种 类 型 ， 受 迫 天 平 动 与 自由 天 平 动 。 受 迫 天 平 动 是 由 于 地 
球 、 太 阳 和 行星 的 引力 引起 的 月 球形 状 上 的 时 变 力矩 ; 自由 天 平 动 主 要 是 地 质 
活动 引起 的 ， 如 冲击 ， 核 - 幅 相 互 作用 ， 或 是 与 受 迫 天 平 动 的 共振 ””。 在 理论 上 ， 
受 迫 天 平 动 的 振幅 、 相 位 和 周期 可 以 根据 月 球形 状 特征 、 弹 性 形变 和 旋转 耗 散 
的 知识 来 计算 。 对 于 自由 天 平 动 ， 只 能 计算 出 周期 ， 必 须 通过 测量 获得 振幅 和 
iu. 

自由 天 平 动 有 三 种 模式 ， 一 种 在 经 度 ， 两 种 在 极 位 。 经 度 模式 是 围绕 月 球 
极 轴 的 周期 为 2. 9 年 的 摆 状 振荡 ”。 剩 下 的 两 种 模式 描述 了 月 球 极 位 与 周期 为 
18. 6 年 的 受 迫 进 动 和 周期 的 
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受 迫 天 平 动 所 决定 的 位 置 之 间 的 变化 。 第 一 个 极 位 模式 被 称 为 摆动 模式 ， 类 似 
于 地 球 的 钱 德 勒 摆动 。 另 外 一 个 极 位 模式 是 自由 进 动 模式 ， 并 对 应 于 空间 中 极 
轴 周 期 为 81 年 的 运动 。 随 着 时 间 的 推移 ， 物 理 天 平 动 会 由 于 能 量 耗 散 而 减弱 
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月 球 天 平 动 可 以 使 用 解析 法 和 数值 法 两 种 计算 方式 ， 在 早期 ， 受 限于 观测 
数据 的 缺乏 以 及 数值 算法 的 计算 量 太 大 ， 数 值 法 相对 于 解析 法 并 没有 太 大 的 优 
势 ， 在 工程 上 对 精度 的 要 求解 析 法 的 摄 动 理 论 完全 可 以 达到 ”。 但 是 随 着 深 空 
探测 的 的 不 断 发 展 ， 对 精度 的 要 求 逐渐 增高 ， 解 析 法 由 于 计算 公式 太 过 复杂 ， 

开始 显示 出 其 应 用 的 劣势 中 。 与 此 同时 ， 在 更 多 的 观测 数据 和 计算 机 技术 的 发 
展 的 支持 下 ， 数 值 法 能 达到 更 高 的 精度 ， 随 着 地 月 数值 历 表 的 建立 ， 使 得 数值 
法 逐渐 成 为 月 球 天 平 动 研究 的 主要 方法 。 月 球 历 表 中 包含 着 月 球 在 国际 天 球 坐 
标 系 中 的 位 置 与 速度 状态 信息 ”"” 。 本 文 对 月 球 轨道 及 物理 天 平 动 的 研究 所 使 
用 的 数据 均 来 自 地 月 历 表 。 

目前 ， 我 国 月 球 探测 已 经 进入 了 月 面 着 陆 的 快速 发 展期 ， 特 别 是 云南 天 文 
台 和 中 山大 学 相继 成 功 开展 了 月 球 激光 测 距 试验 ， 这 为 我 们 加 入 国际 月 球 激光 
测 距 分 析 社 区 提供 了 技术 基础 W”。 随 着 地 面 测 月 站 数量 的 不 断 增 加 ， 国 内 的 
激光 测 月 往 高 精度 的 发 展 ， 未 来 会 积累 大 量 的 高 精度 数据 ， 这 些 数 据 为 历 表 的 
研制 和 物理 天 平 动 的 研究 提供 支撑 。 进 入 新 世纪 后 ， 我 国 提 出 了 月 球 探测 和 深 

宁 测 计划 ， 并 成 功 实施 了 月 球 探测 项 目 。 杨 永 章 等 在 基于 FFT 分 析 的 基础 上 ， 
发 展 了 一 种 频率 提取 和 拟 合 的 数值 方法 ， 从 DE 430 历 表 数据 中 提取 了 月 球 自由 
天 平 动 周期 频率 ， 与 国际 同行 结果 一 致 ， 在 动力 学 模型 方面 ， 建 立 了 基于 国际 
上 最 新 动力 学 模型 的 月 球 自转 双 层 / 单 层 模 型 及 数值 算法 ， 结 果 精 度 与 主流 历 表 
相当 ”。 本 研究 将 实现 使 用 最 新 的 INPOP19a 历 表 实现 对 月 球 物理 天 平 动 的 提取 。 
计算 新 一 代 月 球 历 表 在 物理 天 平 动 数 值 模 型 精度 上 的 提高 ， 为 进一步 高 精度 月 
球 激光 测 距 的 研究 打下 基础 。 


1， 月 球 转动 欧 拉 角 
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月 球 历 表 中 包含 着 各 个 天 体 在 国际 天 球 坐标 系 中 的 位 置 与 速度 状态 信息 ， 
还 包含 了 月 球 的 转动 欧 拉 角 。 利 用 INPOP19a, INPOP17a, DE430 三 个 历 表 中 跨度 
为 600 年 的 数据 提取 得 到 的 欧 拉 角 作 差 比较 如 图 : 
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图 二 INPOP17a 与 DE430 欧 拉 角 的 差别 
BB The difference of Euler angle between INPOP17a and DE430 
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图 三 INPOP19a 与 DE430 欧 拉 角 的 差别 
Fig.3 The difference of Euler angle between INPOP19a and DE430 


结果 如 图 所 示 ，INPOP17 与 DE430 欧 拉 角 差 别 较 大 ，INPOP19 则 更 为 稳定 。 
其 中 对 于 不 同 历 表 之 间 的 差别 稳定 在 0.05 角 秒 来 回 震荡 。 但 是 的 差别 随 着 时 间 
逐渐 增 大 ， 说 明 历 表 中 对 于 的 拟 合 模 型 并 不 稳定 。 对 比 历 表 欧 拉 角 的 差别 ， 计 
算出 地 心 到 月 面 反 射 器 A15 的 距离 最 大 有 30 厘米 的 差别 ， 对 月 球 激光 测 距 的 预 
报 精 度 有 较 大 的 影响 。 对 结果 进行 频谱 分 析 : 
表 1 对 欧 拉 角 进行 频谱 分 析 得 到 的 周期 〈 单 位 : XA) 
Tab.1 The period obtained by performing spectrum analysis on Euler angles (unit: 


day) 


INPOP19a-DE430 INPOP19a-INPOPI17a INPOP17a-DE430 
21.3 27. 32 21.3 
21.3 27. 32 21.3 
21.3 27. 32 21.3 


以 上 过 程 完成 了 对 天 平 动 欧 拉 角 的 提取 ， 并 利用 FFT 算法 对 历 表 中 提取 出 
的 欧 拉 角 之 差 进 行 了 频谱 分 析 ， 获 得 的 结果 如 表 1 所 示 。27.3 天 的 差异 是 极 位 
置 在 惯性 坐标 系 下 自由 天 平 动 的 周期 。 


2， 月 球 天 平 动 与 月 球 转动 的 力学 模型 


月 球 天 平 动 是 对 月 球 运动 在 空间 摆动 的 描述 ， 它 是 月 球 自 转 ， 公 转轨 道 以 
及 周围 天 体 对 月 球 的 力 的 作用 的 综合 表现 ”。 月 球 天 平 动 指 的 是 地 面 观测 者 所 
观测 到 的 月 球 可 见面 上 下 左右 小 幅度 的 摆动 ， 由 于 天 平 动 效应 ， 使 观测 者 能 
到 69% 的 月 面 ， 其 中 19% 时 多 时 人 少 。 

因 月 球 相 对 于 地 心 的 月 面 的 经 度 天 平 动 和 纬度 天 平 动 又 称 光学 天 平 动 或 几 
何 天 平 动 。 几 何 天 平 动 中 的 经 度 项 纬度 项 除了 与 月 球 本 身 的 物理 天 平 动 有 关系 
外 ， 还 主要 与 月 球 的 公转 轨道 有 密切 的 关系 。 因 月 球 实际 摆动 形成 的 天 平 动 则 
称 物理 天 平 动 ， 物 理 天 平 动 是 月 球 自身 相对 于 空间 坐标 的 转动 欧 拉 角 ， 与 月 球 
自身 的 内 部 结构 ， 自 转动 量 ， 转 动 惯量 有 关 。 本 项 研究 具体 实现 从 历 表 中 所 获 
得 的 欧 拉 角 到 月 球 物理 天 平 动 的 转换 过 程 。 

月 固 坐 标 系 分 为 月 球 主轴 坐标 系 〈(Principal Axis Frame, PA) 与 平均 地 球 
指向 坐标 系 (Mean Earth, ME) 。 在 描述 月 球 的 转动 时 ， 我 们 以 月 球 质心 为 原点 ， 
月 球 三 个 惯量 主轴 为 三 个 坐标 轴 ， 建 立 随 月 球 一 起 转动 的 月 球 主轴 参考 系 


(PA) 。 同 时 ， 我 们 可 以 得 到 坐标 指向 与 ICRS 一 致 的 月 心 天 球 参考 系 。 两 个 坐 
标 系 通 过 三 个 欧 拉 角 进 行 转 换 。 
(6) 
这 三 个 欧 拉 角 及 其 变化 率 即 为 月 球 天 平 动 的 数值 表达 ， 反 应 月 球 自 转 的 状 
态 。 依 据 卡 西 尼 定 则 : 月 球 自转 周期 等 于 其 公转 周期 ; 月 球 赤道 与 黄道 的 来 角 
成 ; 月 球 自转 轴 ， 月 球 轨道 平面 的 法 线 以 及 黄道 面 的 法 线 三 者 共 面 ， 且 第 三 者 
处 于 前 两 者 之 间 。 月 球 的 物理 天 平 动 的 理论 基础 ， 主 要 是 刚体 转动 微分 方程 ， 
此 处 方程 由 卡 西 尼 第 一 定律 简化 得 到 ， 对 于 刚体 月 球 转动 微分 方程 为 经 典 的 欧 
拉 - 刘 维尔 方程 : 


在 DE 历 表 中 我 们 可 直接 读 取出 月 心 天 球 参 考 系 到 月 球 主轴 坐标 系 的 转换 矩 
阵 ， 利 用 旋转 矩阵 到 欧 拉 角 的 转换 可 采用 以 下 方法 : 

此 转换 过 程 中 的 旋转 分 别 为 

(8) 
(9) 
(10) 

此 时 有 : 

此 公式 即 为 月 心 天 球 坐 标 系 到 月 球 主 轴 坐 标 系 的 转换 理论 基础 。 

将 这 一 向 量 转换 到 瞬时 黄道 坐标 系 中 : 

(11) 

Hm, B,NP 分 别 代表 坐标 偏差 ， 岁 差 和 章 动 矩阵 。 代 表 瞬 时 真 黄道 坐标 系 
的 旋转 。 

将 PA 坐标 系 中 的 单位 向 量 和 通过 上 述 坐标 转换 过 程 至 瞬时 黄道 坐标 系 中 得 
到 和 。 分 别 代表 在 黄道 坐标 系 下 沿 着 ， 和 的 单位 向 量 。 代 表 坐 标 原点 指向 月 球 
赤道 与 黄道 的 升 交点 的 单位 向 量 。 此 时 有 : 

(12) 


图 四 : PA 坐标 系 到 ICRS 中 的 旋转 欧 拉 角 
Fig.4 Rotated Euler angle from PA coordinate system to ICRS 


图 五 : PA 坐标 系 到 黄道 坐标 系 的 欧 拉 角 


Fig.5 Euler angle from PA coordinate system to ecliptic coordinate system 


利用 图 五 中 的 关系 式 ， 可 以 得 到 : 


(13) 

(14) 

(15) 

(16) 

17) 
根据 Newhall'!: 


(18) 
(19) 
(20) 
在 上 述 过 程 中 ， 我 们 把 图 五 中 赤道 坐标 系 下 的 欧 拉 角 转 换 为 图 六 黄道 坐标 


ARAI. Eckhardt 分 别 使 用 来 表示 天 平 动 中 解 得 周期 项 。 其 中 是 由 分 离 出 来 ， 
被 定义 为 解 的 周期 部 分 ， 由 公式 (190 定义 ， 是 一 个 线性 多 项 式 。 根 据 公 式 
(18) 是 轨道 节点 的 经 度 ， 是 两 个 节点 之 间 的 瞬时 差 ， 将 公式 “18) 变 为 ， 将 
定义 为 以 由 一 个 常数 和 人 解 的 周期 部 分 组 成 。 同 时 “20) 式 可 变化 为 ,定义 由 
一 个 常数 和 结果 的 周期 部 分 组 成 。 至 此 ， 分 别 使 用 来 表示 月 球 的 物理 天 平 动 ， 
根据 上 述 过 程 计算 得 到 结果 最 终 如 图 所 示 : 
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图 六 利用 inpopl9a 历 表 实现 对 月 球 物理 天 平 动 的 提取 
Fig.6 Using inpop19a Ephemeris to Extract the Lunar Physical librations 


至 此 完成 了 利用 数值 历 表 中 欧 拉 角 对 物理 天 平 动 的 转换 。 
利用 上 述 的 方法 ， 从 INPOP17a 5 DE430 中 分 别提 取出 物理 天 平 动 与 INPOP19a 
作 差 比较 结果 如 图 所 示 。 
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图 七 INPOP19a 与 INPOP17a 物理 天 平 动 差别 


Fig.7 the difference between INPOP19a and INPOP17a 
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图 八 INPOP19a 5 DE430 物理 天 平 动 差别 
Fig.8 the difference between INPOP19a and de430 


结果 显示 INPOPI9a 与 INPOP17a 物理 天 平 动 差别 更 为 稳定 。 


3. 物理 天 平 动 的 频谱 分 析 
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由 图 六 物理 量 绘制 的 图 与 图 八 中 INPOP19a 5 INPOP17a 之 差 对 比 发 现存 在 一 个 稳定 频 
周期 ， 且 此 频率 与 本 身 的 固有 频率 不 同 。 故 在 研究 过 程 中 重点 对 此 频率 进行 分 析 。 
首先 ， 根 据 上 述 过 程 ， 利 用 INPOP19a, INPOP17a, DE430 三 个 历 表 中 跨度 为 600 年 的 数 
取得 到 的 物理 天 平 动作 差 比较 如 图 : 
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图 九 INPOP17a 与 DE430 物理 天 平 动 的 差别 
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Fig.9 The difference of physical librations between DE430 and INPOP17a 
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图 十 INPOP19a 与 DE430 物理 天 平 动 的 差别 
Fig.10 The difference of physical librations between INPOP19a and DE430 
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图 十 一 INPOP19a 与 INPOP17a 物理 天 平 动 的 差别 
Fig.11 The difference of physical librations between INPOP17a and INPOP19a 


利用 FFT 算法 对 上 面 的 图 像 进行 频谱 分 析 得 到 如 表 2 的 结果 : 
表 2 对 物理 天 平 动 进行 频谱 分 析 得 到 的 周期 (单位 : XD 


Tab.2 The period obtained by performing spectrum analysis on physical librations 


(unit: day) 
INPOP19a-INPOP17a INPOP19a-DE430 INPOP17a-DE430 
21.21 2T. 18 27.18 
1056. 53 1056. 53 1056. 53 
2T. 21 27.18 21.18 


其 中 1056 天 的 周期 是 经 度 方 向 上 的 自由 天 平 动 周期 ，21. 21 天 和 21. 18 天 则 是 F 和 wobble 


引起 的 。 
4. 总 结 


本 研究 利用 已 有 的 数值 历 表 数 据 ， 开 展 对 月 球 物理 天 平 动 的 研究 。 得 到 


INPOP19a 历 表 中 三 个 欧 拉 角 的 相位 时 间 图 像 ， 与 INPOP17a, DE430 历 表 进行 对 比 ， 
发 现 其 差 值 存在 着 27. 3 天 的 周期 ， 识 别 出 了 极 位 置 在 惯性 系 下 的 自由 天 平 动 。 
进一步 利用 历 表 数据 对 月 球 天 平 动 进行 分 析 ， 根 据 卡 西 尼 定 则 以 及 国际 天 球 参 
考 系 到 月 球 PA 坐标 系 的 转换 模型 ， 提 取得 到 月 球 相对 于 空间 的 摆动 即 月 球 的 物 
理 天 平 动 。 结 果 显 示 ， 月 球 的 物理 天 平 动 存 在 着 一 定 的 周期 ， 但 振幅 在 200 3E 
角 秒 内 ， 与 我 们 能 看 到 的 月 面 是 69% 相 比 ， 是 一 个 极其 微小 的 摆动 量 。 即 在 地 
球 上 观测 月 球 的 天 平 动 时 ， 月 球 的 几何 天 平 动 为 主要 影响 因素 。 把 三 个 历 表 提 
取得 到 的 物理 天 平 动作 差 比较 ， 发 现存 在 一 个 稳定 周期 ， 且 此 周期 与 本 身 的 固 
有 周期 不 同 ， 频 谱 分 析 以 后 得 到 的 此 周期 并 进行 了 比较 。 

对 提取 得 到 的 欧 拉 角 进 行 分 析 ， 计 算得 到 新 一 代数 值 历 表 INPOP19a 与 
DE430，INPOP17a 提取 得 到 的 欧 拉 角 差 别 换算 到 月 球 激光 测 距 上 最 多 有 30 厘米 
的 差别 。 此 差别 对 月 球 激光 测 距 的 预报 精度 产生 较 大 的 影响 。INPOP19a 在 与 
DE430，INPOP17a 对 比 中 ， 显 示 出 较 高 的 稳定 性 ， 故 在 月 球 激光 测 距 或 者 研究 
物理 天 平 动 时 ， 推 荐 采用 历 表 INPOP19a。 


Research on Lunar Physical librations Based on 


Numerical Ephemeris 


Abstract: Lunar librations is a description of the movement of the moon in 
space. To determine the lunar librations makes it possible to speculate about the 
excitation and dissipation mechanisms of the libration, such as meteor strikes, moon 
earthquakes, and viscous friction between the core and mantle. Therefore, the 
determination of the lunar librations is useful for understanding the origin, evolution, 
and structure of celestial bodies in the solar system.During the research, I use the 
latest INPOP2019a ephemeris data, the conversion of Euler angle extracted from the 
ephemeris to the lunar physical librations was completed.The numerical model of the 
physical librations is compared with the ephemeris INPOP17a and DE430.It is found 
that the difference in physical librations of different ephemerides has a stable period 
.It is calculated that there is a maximum difference of 30 cm between the center of the 
earth and the lunar reflector A15. This result has a greater impact on the prediction 
accuracy of the lunar laser ranging.It provides theoretical support for the high- 
precision lunar laser ranging.The results show that the accuracy of INPOP19a is 


significantly higher than INPOP17a and DE430, so INPOP19a is recommended in the 


research of lunar librations. 


Keywords: lunar ephemeris, Euler angle, lunar librations, Fourier transform 
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